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Abstract 



We consider the problem of discretization for U(l) -invariant nonlinear wave equations in any dimension. 
We show that the classical finite-difference scheme used by Strauss and Vazquez 1SV78I conserves the 
positive-definite discrete analog of the energy if the grid ratio is dt/dx < l/y/n, where dt and dx are the 
mesh sizes of the time and space variables and n is the spatial dimension. We also show that if the grid ratio 
is dt/dx = then there is the discrete analog of the charge which is conserved. 

We prove the existence and uniqueness of solutions to the discrete Cauchy problem. We use the en- 
ergy conservation to obtain the a priori bounds for finite energy solutions, thus showing that the Strauss - 
Vazquez finite-difference scheme for the nonlinear Klein-Gordon equation with positive nonlinear term in 
the Hamiltonian is conditionally stable. 

Keywords: Nonlinear wave equation, nonlinear Klein-Gordon equation, discrete space-time, finite-difference 
schemes, grid ratio, U(l)-invariance, energy conservation, charge conservation, a priori estimates. 
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; 1 Introduction 

O, 

We study the U(l)-invariant nonlinear wave equation discretized in space and time. Our objective has been to 
find a stable finite-difference scheme for numerical simulation of the nonlinear wave processes, which corre- 
sponds to a well-posed Cauchy problem and provides us with the a priori energy bounds. 

The discretized models are widely studied in applied mathematics and in theoretical physics. Such models 
originally appeared in the condensed matter theory, due to atoms in a crystal forming a lattice. Now these 
models occupy a prominent place in theoretical physics, in part due to some of these models (such as the Ising 
model) being exactly solvable. Lattice models are also used for the description of polymers. 

The paper [SV78 ] set the ground for considering the energy-conserving difference schemes for the nonlinear 
Klein-Gordon equations and nonlinear wave equations. The importance of having conserved quantities in the 
numerical scheme was illustrated by noticing that instability occurs for a finite-difference scheme which does 
not conserve the energy. The authors gave the implicit difference scheme and wrote down the expression for 
the energy conserved by that scheme. This finite-difference scheme was favorably compared to three other 
schemes in MJV90II . The higher dimensional analog of the Strauss- Vazquez scheme and the corresponding 
energy-momentum tensor was written in [YHHH95]. The general theory of finite-difference schemes for the 
nonlinear Klein-Gordon equation aimed at the energy conservation was developed in [LVQ95] and [FurOl ]. 
In the paper I CJIOl the energy preserving schemes are constructed for a wide class of second order nonlinear 
Hamiltonian systems of wave equations. 
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Importance of the Strauss - Vazquez finite-difference scheme over schemes from [LVQ95 IFurO 1 H is in 
allowing for an easier solution algorithm. Namely, the discrete scheme involves the value of the unknown 
function at the "next" moment of time only at a single lattice point, and can be solved (numerically) with 
respect to the value at that point. See Remark [2721 below. At the same time, the corresponding discrete energy 
for that scheme (see (12.221) below) contains quadratic terms which are not positive-definite, showing that the 
scheme is not unconditionally stable. 



We use the same finite-difference scheme by Strauss and Vazquez MSV78I1 . We show that under the as- 
sumption dt/dx < 1/ y/n on the grid ratio, where n is the number of spatial dimensions, the expression for the 
conserved discrete energy is positive-definite, providing one with the a priori energy estimates in the case of the 
discrete nonlinear Klein-Gordon equation. Moreover, under the assumption dt/dx = 1/ yfn on the grid ratio, 
we show that the equation possesses conserved discrete charge. The continuous limits of our discrete versions 
for the energy and charge coincide with the energy and charge in the continuous case. 

Let us emphasize that the positive definiteness of the energy allows one to have the a priori bounds on the 
norm of the solution. Such a priori bounds are of utmost importance for applications. Numerically, such bounds 
indicate the stability of the finite-difference scheme. 

The discrete charge conservation does not seem to be particularly important on its own, but could be con- 
sidered as an indication that the U(l)-invariance of the continuous equation is in a certain sense compatible 
with the chosen discretization procedure. See the discussion in [LVQ95, Section 1]. 

In Section[2l we formulate our main results. We establish the existence of the solutions to the corresponding 
discrete Cauchy problem and analyze the uniqueness of solutions. We show that for confining polynomial 
potentials the uniqueness will follow when the mesh size is sufficiently small. Besides, we describe a class of 
polynomial nonlinearities for which the value of the mesh size could be readily specified. We also consider the 
charge conservation. The proofs related to the well-posedness are in Section [3j 



2 Main results 
2.1 Continuous case 

Let us first consider the U(l)-invariant nonlinear wave equation 

#(x,t) = A%f>(x,t) -2d x v(x, \^(x,t)\ 2 )i/)(x,t), xeR n , (2.1) 

where ip(x, t) G C^, N > 1, and v(x, A) is such that v G C(R n x R) and v{x, •) G C 2 (R) for each x G R n . 
Equation (12.11 ) can be written in the Hamiltonian form, with the Hamiltonian 



IVf , ivvf 



+ v(x,\i{j(x,t)\ 2 ) 



dx, 



(2.2) 



where for ip G we define \i/j\ 2 = tp ■ ip. The value of the Hamiltonian functional $ and the value of the 
charge functional 

= \ I ■ i> - i ■ il>) dx 



(2.3) 

are formally conserved for solutions to (I2.ll ). A particular case of (12.11 ) is the nonlinear Klein-Gordon equation, 

2 

with v(x, A) = ^A + z(x, A), with m > 0: 



ip = AV> - m 2 i) - 2d x z(x, \ip\ 2 )ip, x G R n , t£ 
In the case z(x, A) > for all x G R n , A > 0, the conservation of the energy J m 



(2.4) 



\4>\ 2 |Vi/>| m \ip 



+ 



+ 



z(x,\4>\ 2 ) 



dx yields an a priori estimate on the norm of the solution: 



(2.5) 
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2.2 Discretized equation 

Let us now describe the discretized equation. Let (X, T) EZ n xZ denote a point of the space-time lattice. 
We will always indicate the temporal dependence by superscripts and the spatial dependence by subscripts. Fix 
e > 0, and let V x {\) = v(eX, A) be a function on Z n x R, so that V x G C 2 (R) for each X G Z n . We 
introduce 

r VxW-VxM A / 

B x {\n) -.= { ' A,/iel, leZ", (2.6) 

[ <9 A Vx(A), A = n, 

and consider the standard implicit finite-difference scheme for (12.11 ) [S V78B : 

^ + '- 2 ff + ^' = £ ^" 2 f + ^" - ifc W + (2.7, 

J'=l 

where ^ 6 C^, JV > 1, is defined on the lattice (X, T) E P x Z. Above, 

ci = (1,0,0,0,... ) G Z n , e 2 = (0,1,0,0,... ) G Z n , etc. (2.8) 

Remark 2.1. The continuous limit of (12.71 ) is given by (12.11 ). with eX corresponding to x G W 1 and tT 
corresponding tot £ R. Since d\V x (\) = B X (X, A), the continuous limit of the last term in the right-hand 
side of (12.71 ) coincides with the right-hand side in (12.11) . 

We assume that ip x takes values in with N > 1. 

Remark 2.2. An advantage of the Strauss- Vazquez finite-difference scheme (12.71 ) over other energy-preserving 
schemes discussed in [ |LVQ95| IFurO 1 II is the fact that at the moment T + 1 the relation (12.71 ) only involves the 
function ip at the point X, allowing for a simple realization of the solution algorithm even in higher dimensional 
case. 

2.3 Well-posedness 

We will denote by ip T the function ip defined on the lattice (X, T) G Z n x Z at the moment T G Z. 
Theorem 2.3 (Existence of solutions). Assume that 

h := inf d x V x (X) > -co. (2.9) 

XeZ™,A>0 

f V-lAi, to < 0; 

T\ = < 

{ +oo, fei > 0. 

Then for any r G (0, ri) a«ci any e > f/iere am/j a global solution ip T , T G Z, ?o Cauchy problem for 
equation (12.71 ) w/?/i arbitrary initial data ip°, ip 1 (which stand for ip T atT = and T = 1). 
Moreover, if^ ,^ 1 ) G Z 2 (Z n ) x Z 2 (Z n ), one has ip T G l 2 (Z n )for allT G Z. 

Remark 2.4. We do not claim in this theorem that ||V ,T |h 2 (z n ) i s uniformly bounded for all T G Z. For the a 
priori estimates on ||V> T ||z2(z™)> see Theorem 12. 1 3 I below. 

One can readily check that any X-independent polynomial potential of the form 

v 

V x (\) = V(\) = Y,C q \ q+ \ C q eR, C p >0 (2.10) 

q=0 

satisfies (12.9b - Note that since liniA^+oo V(X) = +oo, this potential is confining. 
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Remark 2.5. Note that in the case of Vx(X) given by (12.101 ). by the little Bezout theorem, Bx(X, n) defined in 
(12.61 ) is a polynomial of A and \i with real coefficients. 



Theorem 2.6 (Uniqueness of solutions). Assume that the functions iw(A, p) = Bx(X, fi) + 2d\Bx(X, m)(A± 
y/Xfj) are bounded from below: 

k 2 := inf K$(\,n) > -oo. (2.11) 

' \7-lA 2 , &2 < 0; 



T2 



+OO, k 2 > 0. 



Then for any r £ (0, T2) ana* any e > ew'sfs a solution to the Cauchy problem for equation (12.71 ) with 
arbitrary initial data (ip°, ip 1 ), and this solution is unique. 

Remark 2.7. Note that since 

K X (X,X) = B X (X,X) =d x W x (X), 

the values of k\ and k 2 from Theorem 12.31 and Theorem 12.61 whether k 2 > —00, are related by k 2 < fei, and 
then the values of t\ and r 2 from these theorems are related by t 2 <t\. 

Theorem 2.8 (Existence and uniqueness for polynomial nonlinearities). 

(z) The condition (12.111 ) holds for any confining polynomial potential (12. 10b - 

(0) Assume that 

4 

Vx(X)=^2C x , q X q+1 , leP, A>0, (2.12) 

q=0 

where Cx,q > 0/or IgZ™ and 1 < q < 4, ana* Cx,o are uniformly bounded from below: 

k 3 : = inf Cxo > -00. (2.13) 
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fW, ^3 < 0; 

+OO, /C3 > 0. 



Then for any r G (0, T3) ana" any e > there exists a solution to the Cauchy problem for equation (12.71 ) 
with arbitrary initial data (ip ,^ 1 ), and this solution is unique. 

Thus, even though the potential (12.101) satisfies conditions (12.91) and (12.111) in Theorem l2.3l and Theorem l2.6[ 
the corresponding values t\ and t 2 could be hard to specify explicitly. Yet, the second part of Theorem l2.8| gives 
a simple description of a class of X-dependent polynomials Vx(X) for which the range of admissible r > 
can be readily specified. 

We will prove existence and uniqueness results stated in Theorems 12.31 12.61 and !2.8l in Section[3j 
2.4 Energy conservation 

Theorem 2.9 (Energy conservation). Let tp be a solution to equation (12.71) such that i/j t G I 2 {I/ 1 ) for all TeZ. 
Then the discrete energy 



E - (^-?) 2 + Z^4. 4^2 + 2 

xez™ j=l ± 



e n (2.14) 



conserved. 
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Remark 2.10. The discrete energy is positive-definite if the grid ratio satisfies 



r < 1 

e ~ Jn 



(2.15) 



Remark 2.11. The continuous limit of the discrete energy E T defined in (12.141 ) coincides with the energy 
functional (12.21) of the continuous nonlinear wave equation (12.11) . 

Remark 2. 12. If V>° and ij) 1 G Z 2 (Z"), then, by Theorem 1231 one also has ip T G Z 2 (Z n ) for all T G Z as long as 

inf d x V x (X) > -oo. 

xez™,A>o 



Proof. For any u, v € C , there is the identity 



|f | 2 = Re [(u — v ) • (u + u)] 



(2.16) 



Applying (I2.16I ). one has: 



E (l^ +1 -^| 2 -|^-^- 1 | 2 ) = Re E (^ 1 -^- 1 ).(^ 1 -2^+^- 1 )- (2-17) 
Using (12.161 ). we also derive the following identity for any function ip x G C^: 

n 

xez™ i=i 

n 

= Re E E W l ~ S" 1 ) ■ ^x +l ~ 2^ ±ej + ^) + " ^T 1 ) • " Wx +ej + ^) 



Re E 



T+l J.T-1 



1> 



X 



(2.18) 



Further, (12.61 ) together with (12.161) imply that 



(2.19) 



Taking into account (12. 17b . (12.181) . and (12- 19b . we compute: 



xez n 



r 2 e 2 



+ EE 

i=i ± 



Vx(\il> 



X 



4e 2 



Re E 



r+i „7,t-i 



+ 



1 n ^ " 2^ + V'x" 1 



x 



i\t 2 e 2 



n 



+ 



+ 



n(^ +1 - 2^ + i> T x - 1 ) - £ (V^ +e , - 2V£ + ^_ e ,) /T+1 ^ /T _! 



3=1 



+ 



2e 2 2 

The expression in the square brackets adds up to zero due to (12.7b . It follows that E T = E T ~ l for all T G 
Z. □ 
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2.5 A priori estimates 

Theorem 2.13 (A priori estimates). Assume that e > and r > satisfy 



Assume that 



t < 1 
e ~ Jri' 



2 

m 



Vx{\) = ~y X + Wx(X), (2.20) 

where m > 0, and for each I £ Z" the function Wx G C 2 (R) satisfies Wjf(A) > Ofor A > 0. Then any 
solution i[) x to the Cauchy problem (12.71 ) with arbitrary initial data (tp , if; 1 ) G Z 2 (Z n ) x Z 2 (Z n ) satisfies the a 
priori estimate 

H T \\fe n < %, (2.21) 
where E° is the energy (12.141) of the solution ip x at the moment T = 0. 

Proof. This immediately follows from the conservation of the energy (12.141 ) with Vx(A) given by (12.201) . 

pt _ v r/J_ M j^z^g , ^ l^ +1 -^- ej l 2 + l^'-^l 2 

2_> [\ T 2 e 2) 2 ^ 4e 2 

| m 2 (|^ +1 | 2 + |^| 2 ) | iy x (|^ +1 | 2 ) + ^(|^| 2 ) 



□ 

Remark 2. 14. In the continuous limit e — > 0, the relation (12.211 ) is similar to the a priori estimate (12.51) for the 
solutions to the continuous nonlinear Klein-Gordon equation (12.41) . 

Remark 2.15. In HSV78I1 . in the case ij} x G K, (X,T) G Z x Z (in the dimension n = 1), the following 
expression for the discretized energy was introduced: 



o Yl [ T 2 x ^ X+1 — + V{\Vx + \ ) + V{\^x\ )\- (2-22) 



The presence of the second term which is not positive-definite deprives one of the a priori I 2 bound on tp, such as 
the one stated in Theorem l2.13l In view of this, the Strauss- Vazquez finite-difference scheme for the nonlinear 
Klein-Gordon equation is not unconditionally stable. Other schemes (conditionally and unconditionally stable) 
were proposed in [LVO95, [Fur01i Now, due to the a priori bound (12.211) . we deduce that, as the matter of fact, 
the Strauss- Vazquez scheme is stable in n dimensions under the condition that the grid ratio is r/e < 1/ ^fn. 
Note that in the case ifi G M., the Strauss-Vazquez energy (12.221) agrees with the energy defined in (12- 14b . 

2.6 The charge conservation 

Let us consider the charge conservation. We will define the discrete charge under the following assumption: 
Assumption 2.16. 



r __ 1 

e Jn 



(2.23) 



Under Assumption ^. 161 i/j x drops out of equation (12.71) : the latter can be written as 

+ V^' 1 ) (1 + r 2 Bx(\^ x +1 \ 2 , Wx' 1 ] 2 )) = \ E(^ + e, + ^- ej )- (2-24) 



n 
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Theorem 2.17 (Charge conservation). Let Assumption \2J6\ be satisfied. Let tp be a solution to equation (12.241 ) 
such that i)) T £ I 2 {IT 1 ) for allT £ Z (see Theorem \2.3\) . Then the discrete charge 

Q T = i ^ ^ • + • - • ~ • ^-*> n (2 - 25) 

jj conserved. 

Remark 2.18. The continuous limit of the discrete charge Q defined in (12.251 ) coincides with the charge func- 
tional (12.31 ) of the continuous nonlinear wave equation (I2.ll ). 

Proof. Let us prove the charge conservation. One has: 

A n 



xgz™ i=i ± 



4r 



Q T - 1 = E EE ^-cc]=-EEE • ^ - c c. 

xgZ" j=i ± x& n j=i ± 



IS 



n 



Therefore, 

AT{QT ~ n QT ~ 1] = EEE^,r+r)-cc. 

=z™ j=i ± 

E (l + r^d^lM^I^I^ + ^p-CC. =0. 

XGZ' 1 

To get to the second line, we used the complex conjugate of (12.24b - This finishes the proof of Theorem l2.17l □ 

3 Proof of well-posedness and uniqueness results 

Proof of Theorem 1231 We rewrite equation (12.71 ) in the following form: 

9 n 

T 



3=1 

(3.1) 

By (12.91 ) and the choice of ri in Theorem 12. 3 1 for r G (0, ri) one has 

inf (l + r 2 a A Fx(A)) >0. (3.2) 

Since 

inf B x (\,fi)= inf V^(A)-Vy(/x) = ^ q y {x) (3 3) 

XgZ™,A>0,M>0 XeZ",A>0,/j>0,A^/i A — /X XeZ™,A>0 

inequality (13.21 ) yields 

c:= inf (1 + T 2 B x {\fi)) > 0. (3.4) 

XeZ",A>0,At>0 v 7 

Let us show that equation (13.11 ) allows us to find ipx +X , for any given IgZ" and T £ Z, once one knows 



V> T and ^ T 1 . Equation (13.11 ) implies that 

(1 + T 2 B x {\Vx +l \\ \V x - l \ 2 )) + flT 1 ) = fx, (3-5) 



o n 



tx ■= "J E (^+e. " 2 ^ + ^-e,) + 2^ E C*. (3.6) 
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If £ x = 0, then there is a solution to (13.51 ) given by ip x 



l,T+l 



ip x 1 . Due to (13 -4b . this solution is unique. Now 

(3.7) 



let us assume that £ x ^ 0. We see from (13.51) that we are to have 

mil rp -i rp 

4> x + ip x ~ = s£ x , with some s G 
Let us introduce the function 

f(s) := (1 + T 2 B x (\s& ~ 4> T x - l \\ H T x ~ l \ 2 ))s. (3.8) 

T~i 1 rp rp 

We do not indicate dependence of / on tp x , £ x , and X, treating them as parameters. For t^ x ^ 0, we can 
solve (13.51 ) if we can find seR such that 

f(s) = 1. (3.9) 

Since /(0) = 0, while lim^oo f(s) = +oo by (13.41 ). one concludes that there is at least one solution s > to 
(EU). 

Let us prove that once (tp , 6 Z 2 (Z n ) x / 2 (Z n ), then one also knows that ||V ,T |b 2 (z n ) remains finite (but 
not necessarily uniformly bounded) for all T 6 Z. As it follows from (13.41 ) and (I3.5I ). 



|V£ +1 |<^xl + l^ 



1 



T-li 



(3.10) 



Since \\(? 



IP(Z») < ( + 2 



U 2 (Z«) 



by (13.6I ). the relation (13.101 ) implies the estimate 



i,T+1i 



U 2 (Z») 



< 



1 At 2 



c \ e 



+ 2 



u 2 (z«) + 



i.T-li 



u 2 (Z")> 



(3.11) 



and, by recursion, the finiteness of ||^ T ||/ 2 (z™) f° r ai l T > 0. The case T < is finished in the same way. □ 
Now we turn to the uniqueness of solutions to the Cauchy problem for equation (12.71) . 



Proof of Theorem 1261 First, note that, by Remark 12771 t\ from Theorem 12.31 and T2 from Theorem 12.61 are 
related by < t\. Therefore, the existence of a solution ip x to tne Cauchy problem for equation (12.71) follows 
from Theorem 12.3 1 

Let us prove that this solution ip x is unique. When in (13.61) one has 



(T ._ T _ 



then, by (13.41) . the only solution tpj^ 1 to (13.51 ) is given by ip x +L = —ip x ~ L . We now consider the case £ x ^ 0. 
By (I3.5I ). (I3.7I ). and (I3.8I ). it suffices to prove the uniqueness of the solution to (13 -9b - This will follow if we show 
that f(s) satisfies 

f'(s)>0, seR. (3.12) 

The explicit expression for f'(s) is 

1 + t 2 2?x(|s£ " lift- 1 ! 2 ) + ^A^xd^ - ^r 1 ! 2 , l^" 1 | 2 )(-2Re • £) + m T x\ 2 s)s. 

(3.13) 

Using the relation (13.7I ). we derive the identity 



2 ™ 



£W£+e, " 2 ^X + *rX- ej ) + 2^5 = 0, 
T+l _ „/.T-l 



(-2Re(^~ 1 -^) + 2|ei|^)s = 2 
and rewrite the expression (13.131) for /'(s) as 



T 



2 l^" 1 ! 2 



th T - 1 



+ 



1^ 



T-l 1 2 
X 



/'W = l + r 2 [i?x(|^ +1 | 2 ,l^- 1 | 2 ) + 25 A ^(|^ +1 | 2 ,|^- 1 | 2 )(l^ 



x + 2 I - 



T-l |2, 
X 



(3.14) 
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We denote A = |V>x +1 | 2 , P = IV'x" 1 ) 2 - Since A - y/Xp, + f < + ^^H 2 ^ A + + 4> we see that 

f'(s) > 1 +r 2 min^ ^inf > if|(A,/i), with K|(A, /j) = 5 X (A, /j) + 2d x B x {X, /i)(A ± v 7 ^)- 

~" (3-15) 
By (12.1 II) and by our choice of T2 in Theorem 12. 61 for any r G (0, T2) we have 

k:= inf {l + t 2 K±(\,u)\ > 0; 

then, by (13.151 ). /'(s) > ft, where k > 0. It follows that for (| / there is a unique s which solves (13.91) . 
Hence, there is a unique solution to e q uat ion (13.51 ) for given values ip x 1 an d £x- This finishes the proof 
of the Theorem. □ 

Proof of Theorem [2751 Let us prove that the condition (12.111) in Theorem l2.6l is satisfied by any polynomial po- 
tential of the form (12.101) . The inequality (12.1 II ) will be satisfied if the highest order term from V(X) contributes 
a strictly positive expression. More precisely, we need to prove the following result. 

Lemma 3.1. Let V(X) = X p+1 , so that B(X,p) = AP ^Zjf — , p > 0. Then the following inequality takes 
place: 



inf 

A>0,a*>0, A 2 +^i 2 =l 



B(X, p) + 2d x B(X, p)(X ± VA/i) 



>0. (3.16) 



Proof. Since B and d\B are strictly positive for A 2 + /i 2 > 0, the inequality (13.161 ) is nontrivial only for the 
negative sign in (13.161 ) and only when p > A. First we note that 

r-% / \ \ ii p+1 -\ p+l a of . , -(p+l)\P(v-\)-Xf> +1 + ^ +1 pP +1 -(p+l)\Pp + p\P +1 
B{X ^ )= p-X ' 9 ^^ = ^a)2 = {^rxf 

Let z > be such that z 2 = X/p. To prove the lemma, we need to check that 



1 - z 2 p+ 2 Hp + 1)^ + P z 2 p+ 2 ( 2 
1-z 2 + (l-z^ 2 



+ 2 ^ ' ' ' ^ {z z -z)>0, 0<z<l, (3.17) 

L — Z~ 

or equivalently, 



(1 + z)(l- z 2p+2 ) - 2z(l - (p + l)z 2p + pz 2p+2 ) > 0. 
The left-hand side takes the form 

(1 + z )(l - z 2 p+ 2 ) - 2z(l - z 2p+2 - (p+ l)(z 2p - z 2p+2 )) = (1 - z){\ - z 2p+2 ) + 2z(p + l)(z 2 P - z 2p+2 ), 
which is clearly strictly positive for all < z < 1 and p > 0, proving (I3.17I ). □ 

This finishes the proof of the first part of Theorem 12. 8 1 now we turn to the second part. 

Lemma 3.2 (Uniqueness criterion). Assume that for a particular r > and for all X > 0, p > 0, X £ IT 1 , the 
following inequalities hold: 

1 + r 2 inf (B x (X,p)-d x B x (X,p)^) > 0; (3.18) 
inf 5 AJ B x (A,/i) >0. (3.19) 

XeZ™, A>0,ai>0 

77je?i jj a solution -ip x to the Cauchy problem for equation (12.71) with arbitrary initial data (t/j ,^ 1 ), and 
this solution is unique. 
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Proof of Lemma [Ol The inequalities (13.181 ) and (13.191 ) lead to 



1 + r 2 inf B x {X,X)>0, 

X£Z n ,\>0 

hence, by the same argument as in Theorem l2.3l there is a solution ip x . The relation (13.141 ) shows that f'(s) > c 
for some c > 0. The rest of the proof is the same as for Theorem 12.61 □ 

In the second part of Theorem 12. 8 1 we assume that 

4 

VxW = J2Cx, q X q+ \ ia n , A>0, (3.20) 

q=0 

where C x ,q > for X G Z n and 1 < q < 4, and 

k 3 = inf Cxo > -oo. (3.21) 

One can see that the term Cx,oX in Vx(X) contributes to Bx(X,fi) the expression bx,o(X, fi) = Cx,o, while 
each term in Vx(A) of the form Cx,qX q+1 , with 1 < q < 4 and Cx,q > 0, contributes to Bx(X,fj,) the 
expression Cx, q b q (X, n), with 6 g (A, = ^fe=o X q ~ k /J, k . For r G (0, T3), with T3 = yj — I/&3 for £3 < and 
T3 = +00 for &3 > 0, one has 

1 + r 2 inf C x > 0. (3.22) 

Lemma 3.3. For 1 < q < 4, 6 g (A, //) = Xlfc=o X q ~ k /J, h satisfies the inequality 

b q (X,fi) > d\b q (\,n)^ for all A, fi > 0. 

By (13.221 ) and Lemma 13.31 condition (13.181 ) is satisfied. Since Cx, q > for 1 < q < 4, each term 
Cx,qb q (X, /i) satisfies condition ( 13.191) . Therefore, by Lemma [3721 there is a unique solution t/j x to the Cauchy 
problem for equation (12.7b - This finishes the proof of Theorem 12. 81 □ 



4 Conclusion 

We found out that the Strauss - Vazquez finite-difference scheme HSV781I for the U(l)-invariant nonlinear wave 
equation in n spatial dimensions, with the grid ratio r/e < 1/ yfn, admits the positive-definite discrete analog 
of the energy which is conserved. The result holds in any spatial dimension n > 1, for the field valued in 
C N , N > 1. In the case of the nonlinear Klein-Gordon equation with positive potential, this provides a priori 
bounds for the solution, showing that the finite-difference scheme is stable. 

We found out that if the grid ratio is r/e = 1/ y/n, then this finite-difference scheme also preserves the 
discrete charge. 

We proved that the solution of the corresponding Cauchy problem exists and is unique for a broad class 
of nonlinearities. In particular, this is the case for any confining polynomial potential if the discretization is 
sufficiently small. Finally, we indicated a class of polynomials for which the size of the discretization could be 
readily specified. 
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